Low-energy sector of the S=l/2 Kagome antiferromagnet 
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Starting from a modified version of the the S=l/2 Kagome antiferromagnet to emphasize the role 
of elementary triangles, an effective Hamiltonian involving spin and chirality variables is derived. 
A mean-field decoupling that retains the quantum nature of these variables is shown to yield a 
Hamiltonian that can be solved exactly, leading to the following predictions: i) The number of 
low lying singlet states increase with the number of sites A'^ like 1.15^; ii) A singlet-triplet gap 
remains in the thermodynamic limit; iii) Spinons form boundstates with a small binding energy. 
By comparing these properties with those of the regular Kagome lattice as revealed by numerical 
experiments, we argue that this description captures the essential low energy physics of that model. 
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Despite a very intense activity over the past 10 years, 
the magnetic properties of the S=l/2 Kagome antiferro- 
magnet remain an open problem. If a number of facts 
seem to be rather firmly established by now thanks to 
the very extensive numerical simulations that have been 
performed on that system [0-^ , a simple theoretical pic- 
ture that accounts for the basic findings has not emerged 
yet. The most striking feature is probably the presence 
of many, low-lying singlet states The first indica- 
tion that this might be the case was the appearance of 
a low temperature peak in the specific heat. While the 
evolution of this peak with the size of the system is not 
clear yet, the numerical determination of all the low-lying 
singlet states for systems with up to 36 sites shows that 
their number increases like 1.15^, where N is the num- 
ber of sites of the system. The best candidate to explain 
this proliferation of low-lying singlets is a short-range 
RVB description of the low-energy sector based on dimer 
coverings of the Kagome lattice with nearest-neighbour 
singlets |9- 
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The main problem with this approach 
is that the number of dimer states increases like 1.26^, 
i.e. much too fast ]l^, and no convincing criterion could 
be found that allows one to select the relevant singlet 
states. The other important, although less accurately es- 
tablished, findings of the numerical simulations are the 
absence of long-range magnetic order in the ground-state 
and the presence of a singlet-triplet gap in the thermody- 
namic limit 1^. Finally the role of spin 1/2 excitations, as 
well as the consistency of the numerical results with some 
exotic types of order [^3]-^ , is still under investigation. 

In this paper, we propose a simple explanation of these 
properties. We start from the following observation: The 
exponential increase of the number of these low-lying 
states suggest that they originate from the partial lifting 
of a local degeneracy that would be present if some of 
the couplings were set to zero. Now the natural bricks to 
construct the Kagome lattice are triangles, and spins 1/2 



on a triangle lead to a fourfold degenerate groundstate: 
two doublets that differ by their chirality. So let us inves- 
tigate how this degeneracy is lifted if one constructs the 
Kagome lattice by coupling triangles. This amounts to 
studying the modified Kagome lattice depicted in figure 
1 starting from the limit J'/ J ^ 1. This can be seen 
as a triangular lattice of triangles with Nt — N/3 sites, 
where is the number of sites of the Kagome lattice. 



FIG. 1. Sketch of the dimerized Kagome lattice, 
lines: J, dashed lines: J'. 
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The first step is to derive an effective Hamiltonian 
in the subspace of the groundstates of the triangles, as 
in Subrahmanyam's block spin perturbation approach 
to the non-dimcrized Kagome lattice jl^]. Follow- 
ing Schulz's approach to the problem of three coupled 
Heisenberg chains with periodic boundary conditions 
jlSf , we describe the four groundstates of a triangle with 
two Pauli matrices: a for the spin of the doublet, the 
eigenstates of CTz being denoted f and | , and r for its 
chirality, the eigenstates of being denoted R for right 
and L for left. In terms of the original spins S, these 
states can be written 
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\aR >— —7={\ — aaa > +uj\a — aa > +uj'^\aa — a >) 
v3 

\aL >— —^{\ — aaa > +uj'^\a — aa > +uj\aa — a >) (1) 
v3 

where u = exp(27ri/3) and a =| or I. \a1a2a3 > rep- 
resents a configuration of the original spins S within one 
triangle, the indices corresponding to the convention of 
Fig. 1. Note that the total spin is now given by ■ (7^)^, 
where the prime means that the sum runs over the tri- 
angular lattice. Each triangle has an energy — 3J/4, and 
energies will be measured with respect to the groundstate 
energy -(3J/4)iVt of the J' /J = case. Then, to first 
order in J' , the effective Hamiltonian H on the triangular 
lattice is given by: 

H = (J79) H^HJ^, ^ - 



<4J> 



<4J> 



Hl^ = (1 - 2(a,,r,r + 4t+))(1 - 2(A,rr + ^r/)) (2) 

where < i,j > denotes pairs of nearest neighbors. In 
H^j , aij and Pij are complex parameters that depend on 
the type of bond: (resp. Pij) equals 1, or a; when 
the original spin in triangle i (resp. j) involved in the 
bond sits at site 1, 2 or 3 with the convention of 

Fig. 1. In the basis \RR >, \RL >, \LR > and \LL > 
the eigenstates of H^j can be easily calculated: 

£'1=9 

\<p2{h3) >= ^{^,f3ij,aij,a^jf3,j) E2 ^ I 
I'l^Uhj) >= = -3 

\(|}l{hj)>^^{'^,P^j,-a^J,-a^jPiJ) E4 = -3 (3) 

while the eigenstates of H[j are denoted \4>'s,m{'''' j) '^i^h 
energies —3/4 for the singlet {S = 0, m = 0) and 1/4 for 
the triplets (5 = 1, m = 0, ±1). 




FIG. 2. Solution of the two-triangle problem. Singlets are 
represented as double lines. 

It is useful to start with two triangles coupled by a sin- 
gle link. In the spirit of the Majumdar-Ghosh solution of 
the Ji — J2 chain [|^, the groundstate wavefunction in 
terms of the spins 5* is obtained as the product of three 
singlets involving respectively the link between the tri- 
angles and the remaining two spins on each triangle (see 
Fig. 2), and its energy is — 6J/4— 3J'/4. In that simple 
case, the effective Hamiltonian Hij is the product of a 



spin part H[j and a chirality part H^j, and its eigenval- 
ues are the products of the eigenvalues of H^j with those 

of H(y The groundstate energy of Hij is thus — 3J'/4: 
The exact groundstate energy is recovered, as it should 
since it depends linearly on J'. 

In the general case of Eq.(|^), the Hamiltonian is not 
the product of a spin part and a chirality part, and its 
solution is in principle as difhcult as the original prob- 
lem. However, we note that the Hamiltonian of Eq.(||) is 
formally similar to the Kugel-Khomskii model that was 
introduced in the context of orbitally degenerate magnets 
p2| , and a mean- field decoupling of the spin and orbital 
degrees of freedom is known to give an accurate picture 
of the physics when the assymmetry between spin and 
orbital degrees of freedom is strong enough psf . In our 
case, the assymmetry between spin and chiral variables 
in Eq.(|^) is clearly very strong, and such a decoupling 
amounts to the replacement of H with the mean-field 
Hamiltonian: 



H 



MF 



E 



(4) 



where the parameters ajj =< HJj > and afj =< Hf^ > 
have to be determined self-consistently. Note that this 
mean-field problem is still very complicated a priori since 
it involves S=l/2 Heisenberg like models on a triangular 
lattice. 

Remarkably enough, the low-energy solutions of that 
problem can be determined analytically. Let us concen- 
trate for the moment on clusters with an even number 
of sites and with periodic boundary conditions, and let 
us consider a dimer covering of the triangular lattice 
by nearest-neighbor dimers. Denoting by D the set of 
nearest-neighbor pairs that enter this covering, we can 
construct a wavefunction |<I'o(^) > in the following way: 

\MD)>^ n Wi{hj)>mifl{hj)> (5) 

<jj>e-D 

Clearly |$o(-D) > will be a solution of the problem if 
alj = afj = as soon as < i,j >^ D. This turns out to 
be true thanks to the following properties 

<q^l{i,j)ct>l{k,l)\Hj^\ct>l{i,m{k,l) > =0 
< ro,oihjWo,oikJ)\HJk\ro,oihjWoA>'J) > - (6) 

which can be easily checked directly with the expressions 
of the wavefunctions. So \^o{D) > is a solution charac- 
terized by ajj = 9, afj = -3/4 if < i,j >e D and oth- 
erwise, its energy is given by Eo{S = 0) — (3J'/8)iVt, 
and it is a singlet p4| . To prove that is minimizes the 
energy of Eq.(^, we have solved the mean- field prob- 
lem numerically on small clusters with up to 3 x 4 sites. 
This has been done by iteration starting from random 

values of of; and , and we found that the lowest en- 

*j ''J 

ergy is always equal to — (3J'/8)A^t. In fact, solutions like 
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^o{D) exist most of the time for such mean-field Hamil- 
tonians, but they are usually significantly higher in en- 
ergy than uniform solutions. In the present case however, 
other solutions involving \ (f>2{h j) > sre very bad energet- 
ically because E2 — Ei/9. Note that the wavefunction 
> corresponds neither to ferromagnetic nor an- 
tiferromagnetic ordering of the chiral variables since it 
is a linear combination of the four basis states \RR >, 
\RL >, \LR > and \LL >. So there is no chiral ordering 
of the type discussed by Baskaran for the triangular 
lattice even locally. 

Similarly, triplet solutions can be constructed for a 
given dimer covering of the triangular lattice. Choose 
two neighbouring sites (ioi jo)j and consider a dimer cov- 
ering D{iQ,jo) of the remaining sites. The wavefunction 
with lowest energy is of the form 

\MD{ioJo)) > (IC(*o,Jo) > <»|0M„(«o, jo) >) (7) 

In this expression, m — 0, ±1 and n can take the values 3 
or 4 H. The energy of this state is Eo{S = 1) = Ea{S = 
0) -I- (2/3) J'. Again it was checked numerically on a 3 x 4 
cluster that this is indeed the lowest energy in the triplet 
sector. So this mean- field approach predicts that there is 
a singlet-triplet energy gap A equal to (2/3) J'. 

Another class of low-lying states exist in the triplet 
sector. They can be constructed in the following way: 
Choose two sites {k,l) that are not nearest neighbours, 
denote by D(k, I) a dimer covering of the remaining sites, 
and consider the wavefunction 

\^o{D{k,l))>\<JkTk>\aiTi> (8) 

where \akTk > (resp. \aiTi >) can be any configuration at 
site k (resp. I). Then similar arguments show that this is 
a solution with energy Ei{S = 1) = Eo{S = 1) + J'/ 12. 
Each unpaired site correspond to a S' = 1/2 excitation 
and can be interpreted as a spinon. This mean-field ap- 
proach predicts that spinons form triplet bound states on 
neighbouring sites with a binding energy equal to J'/12. 
In fact, if we consider a cluster with an odd number of 
sites, the groundstate can be shown to consist of one un- 
paired site - i.e. one spinon - times |$o(-D) >, where D 
is a dimer covering of the remaining sites. 

Now let us turn to the very interesting question of the 
groundstate degeneracy. The energy Eo{S = 0) does not 
depend on the particular dimer covering D used to con- 
struct the wavefunction |<I>o(I?) >. So, for a given cluster, 
the degeneracy is controlled by the number of dimer cov- 
erings. This number can be calculated using standard 
techniques |2^,^. For the triangular lattice, we found 
that it increases with the number of sites Nt like a^* with 
\na — lo^ lo^ \n{4: + 4 sin x sin y + 4: sin^ y)dxdy. A 
numerical integration yields In a = 0.4286, or a = 1.5351. 
In terms of the original Kagome lattice, this corresponds 
to a degeneracy that increases like {a^/^)'^ = 1.1536^ 
since N = SNt. 

Next let us try to assess the validity of this mean-field 
approach. One argument in favour of this kind of ap- 
proach is that it is qualitatively correct in the case of 



three coupled chains with periodic boundary conditions: 
Schulz's analysis based on Renormalization Group argu- 
ment predicts a dimerized, twofold degenerate ground- 
state with gapped spin excitations ||l^ , |20|] , and this is 
exactly the physics of the mean-field solution adapted to 
that case To check its validity in cases of higher 

groundstate degeneracy, we have performed numerical 
simulations of the effective Hamiltonian of Eq.(H) for 
finite systems of coupled triangles. For the modified 
Kagome problem this turns out to be very difficult be- 
cause the first cluster where dimer coverings can be done 
properly has 3x4 sites. To determine the low lying states 
is then a numerical task comparable to the Kagome prob- 
lem itself for 36 sites, and this is beyond the scope of the 
present paper. So we have compared mean-field and ex- 
act results for similar but different systems, namely for 
rings of n triangles with diagonal bonds up to n — 6. 
For n = 2, the mean-field description is of course ex- 
act. For n — 4 and 6, the structure of the low-energy 
spectrum is correctly given by the mean-field approach. 
For instance, for n = 4, the mean-field approach pre- 
dicts that the groundstate is threefold degenerate, and 
the exact results give 3 low-lying singlets at — 1.72053J' 
(non degenerate) and —1.54077 J' (two-fold degenerate), 
the next state being a triplet located at —1.13942 J'. For 
n = 6, the groundstate degeneracy is 6 at the mean-field 
level, and the exact results also give 6 low- lying singlets. 
The details will be given elsewhere, but the answer is 
clear: The mean-field approach is qualitatively correct 
in predicting the number of low-lying singlets. To go 
beyond mean-field is expected to partially lift the degen- 
eracy within the groundstate manifold but not to change 
the number of low- lying singlets. 

Finally, let us come to the most important question: 
What can we learn from this approach concerning the 
regular (non-dimerized) Kagome lattice? First of all, 
let us translate our results in terms of the original spins 
S. We know from the analysis of the two triangle prob- 
lem that the basic brick of our mean-field wavefunction, 
namely a two-site wavefunction of the type \(f>l(i,j) > 
(8)|0Q g(z, j) >, corresponds to a dimer mapping of the two 
triangle problem (see Figure 2). So, the wavefunctions of 
Eq.(^ correspond to a certain subset of the dimerized 
wavefunctions of the Kagome lattice used by Zeng and 
Elser |ll). Next, let us compare the physical proper- 
ties of this mean-field solution with the numerical results 
obtained on the Kagome lattice. As we noticed in the 
introduction, the number of dimerized states increases 
like 1.26^, i.e much too fast. However, the number of 
states selected with our criterion increases roughly like 
1.15^, in agreement with the numerical results for even 
clusters This is probably the most interesting result 
of the present approach: It provides a simple but nev- 
ertheless quantitative explanation of the very numerous 
low- lying singlet states of the S=l/2 Kagome antiferro- 
magnet. Their number scales with the number of dimer 
coverings of the underlying triangular lattice of triangles. 
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The fact that one can choose triangles pointing upwards 
or downwards to build the wave-function means that the 
actual low lying states can have domains with different 
orientations of the triangles. Elementary energy consid- 
erations suggest that these domains are large however, 
so that they can only marginally contribute to the in- 
crease of the number of singlet states with the size of the 
system. Interestingly enough we can also explain the ap- 
parent discrepancy between odd and even clusters in the 
results of Ref. jsj. For odd clusters, there is an unpaired 
site, and according to the present theory, the degeneracy 
is expected to scale like N x 1.15^. Omitting this pref- 
actor in the fit leads to an overestimate of a (see Figure 
3). 



Zurich where this project started. The numerical simula- 
tions were performed on the Cray supercomputers of the 
IDRIS (Orsay, France). 



FIG. 3. Logarithm of the number of dimer coverings of the 
triangular lattice as a function of the number of sites Nt for 
small systems with periodic boundary conditions. The slope 
of \n{Nodd) is 1.18, while the slope of both \n{Neven) and 
HNodd/Nt) is 1.15. 

To summarize, we have studied a dimerized version 
of the S=l/2 Kagome antiferromagnet by performing a 
mean-field analysis of the effective hamiltonian that de- 
scribes its low-energy sector. This approach leads to a 
transparent picture of the low energy properties which 
turns out to bear remarkable similarities to those re- 
ported for the regular Kagome lattice e.g. the number 
of low-lying singlets or the presence of a singlet-triplet 
gap. Besides, these results lead to natural subspaces of 
dimer wavefunctions to describe the low energy singlet 
and triplet sectors of the regular Kagome model. A vari- 
ational study of this model using these wave-functions is 
in progress to make more precise statements about the 
accuracy of this description of the low energy sector of 
the regular Kagome antiferromagnet. It is expected to 
give useful information concerning the structure of the 
low- lying singlet sector, for instance the low temperature 
specific heat, and to allow one to make more precise pre- 
dictions concerning the persistence of the singlet-triplet 
gap and the binding of spinous in the thermodynamic 
limit. 
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